The singular field used to calculate the self-force on non-spinning and spinning 

particles 
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The singular field of a point charge has recently been described in terms of a new Green's function 
of curved spacetime. This singular field plays an important role in the calculation of the self-force 
acting upon the particle. We provide a method for calculating the singular field and a catalog of 
expansions of the singular field associated with the geodesic motion of monopole and dipole sources 
for scalar, electromagnetic and gravitational fields. These results can be used, for example, to 
calculate the effects of the self-force acting on a particle as it moves through spacetime. 



I. INTRODUCTION 

Detailed knowledge of the evolution of a binary sys- 
tem consisting of a stationary supermassive black hole 
(of mass of the order of 1O 6 M0) and a stellar-mass neu- 
tron star or black hole would be an important aid for 
the analysis of data from space-based gravitational wave 
detectors such as LISA In general, the stellar- mass 
neutron star or black hole, usually modelled as a point- 
like particle, inspirals toward the central black hole due 
to radiation reaction. Radiation reaction is only one con- 
sequence of the self-force on a particle, namely the force 
which results from the interaction of the particle with its 
own gravitational field. 

Generally, the spin of the stellar-mass compact object 
is expected to have an effect on the gravitational wave- 
forms by such a system. That effect is expected to be 
small and does not need to be taken into account in 
template waveforms aimed at detection of such binaries. 
However it will have to be taken into account for accurate 
estimation of the physical parameters of the binary sys- 
tem. For that reason, self-force calculations for combined 
monopole and dipole sources are required. This paper is 
a step toward self-force calculations for spinning sources. 



A. Self-Force Calculation 

One method for calculating the self-force on a parti- 
cle generating an electromagnetic field was suggested by 
DeWitt and Brehme That method involves decom- 
posing the retarded field at point p, generated by the 
charged particle at point p', into its direct part (which 
comes from the part of the Green's function with support 
only on the null cone of point p) and its tail part (which 
comes from the part of the Green's function with support 
only within the null cone of point p) , and then using the 
tail part of the electromagnetic field in the Lorentz force 
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law. Similar analysis by Mino, Sasaki and Tanaka [|[ 
described the gravitational self-force on a massive par- 
ticle in terms of the tail part of the gravitational field. 
Axiomatic approaches were used by Quinn for the scalar 
field [i| and by Quinn and Wald for the electromagnetic 
and gravitational fields [f| and arrived at similar conclu- 
sions. 

The actual calculation of the self-force on a particle 
in geodesic motion around a Schwarzschild black hole, in 
terms of the tail part of the field, has been successfully 
performed by a variety of groups [g, 0, H, H, EH EH EH ES 
Il4 Il5j . In these applications, the retarded field is first 
calculated via separation of variables coupled with a sum 
over modes. Analytic methods are used to calculate the 
part of the self- force from each mode of the direct part of 
the field. Then, for each mode, the part of the force from 
the direct part is subtracted from the part of the force 
from the retarded part. Finally, this difference is summed 
over all modes and this sum provides the self-force which 
results from the tail part of the field. It is clear, then, 
that the tail part provides a useful calculational avenue 
to the self-force. However, a few shortcomings remain 
on a fundamental level [if}. The Green's function of 
the tail part has support inside the past null cone of the 
field point, which implies that this description of the self- 
force formally depends upon the entire past history of 
the particle. Also, in some circumstances the tail part is 
finite but not differentiable at the location of the particle 
and averaging of the derivative over different directions 
of approach to the particle is required to determine the 
actual self-force. Perhaps most critically, the tail part of 
the field is not a solution of the field equations. Thus, 
the tail part of the field of a particle can be described 
mathematically, but it cannot by itself be considered a 
field. 



1. Self-force calculation using the singular field 

A different method for calculating the self-force was 
proposed by Detweiler and Whiting [16] for the scalar, 
electromagnetic and gravitational fields. That method 
uses a new Green's function for curved spacetime to de- 
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termine the singular field of the particle. The singular 
field was first defined in [16J where it was shown by con- 
struction to be the part of the retarded field of the parti- 
cle that exerts no force on the particle itself. The prop- 
erties of the singular field were detailed in section 13.5 of 
[l7l ]. For completeness, we give here the definition and 
properties of the singular field. For a more detailed dis- 
cussion, including the motivation for the definition of the 
singular field, we point the reader to [l6| and [l7| . 

Assume a particle of scalar charge q at point p' in 
a background spacetime. We want to calculate the re- 
tarded field of that charge at point p in spacetime. The 
retarded Green's function G ret obeys the inhomogeneous 
equation 

V 2 G* ct (p,p') = -4ttS(p,p') (1) 

where V denotes differentiation with respect to the back- 
ground metric. The retarded Green's function generates 
a field which is the physically relevant solution, because 
it obeys appropriate outgoing-wave boundary conditions 
at infinity. 

The singular Green's function G s is defined [H, [l7| 
as the function that satisfies the same inhomogeneous 
equation as the retarded Green's function 

V 2 G s (p,p') = -4wS(p,p'), (2) 

with the additional requirements that it is symmetric in 
its arguments p and p' and that it vanishes if p is within 
the past or future null cones of p' . By definition, the 
singular Green's function has the same singular behavior 
at the location of the particle as the retarded Green's 
function: the divergence of the field near the particle 
is insensitive to the boundary conditions at infinity. It 
has been shown that such a function always exists in the 
local convex neighborhood of p' (see chapter 13 of [l7j ) 
and that the singular field that results from it does not 
exert any force on the particle itself. As explained in 
[T7I ] . because the singular Green's function is symmetric 
in its arguments, it does not distinguish between past 
and future and it produces a field that contains equal 
amounts of outgoing and incoming radiation at infinity. 
Removing G s from the retarded Green's function G rct 
has the effect of removing the singular behavior of the 
field without that affecting the motion of the particle. 

The radiative Green's function, also referred to as the 
regular remainder, is defined by 

G R (p, P ') = G Ict (p,p')-G s (p,p'). (3) 

This definition implies that the radiative Green's func- 
tion obeys the homogeneous equation 

V 2 G R (p,p')=0, (4) 

that it is equal to the retarded Green's function if p is 
within the future null cone of jpf and that it vanishes if p 
is in the past null cone of p' [17| ■ The radiative Green's 
function produces a field that is smooth and differentiable 



on the world line of the particle and which properly en- 
codes the outgoing-wave boundary conditions at infinity. 

As has already been mentioned, after the singular field 
is subtracted from the retarded field, the radiative field 
is entirely responsible for the self-force. Specifically for 
the scalar case 

TP R (p)=^ et (p)-i> S (l>) (5) 

and the self-force is calculated by using the radiative field 
in the self-force equation 

T a = q lim V a i! K {p) = q lim V a [^ ct {p) - ^ s (p)}. (6) 

p^p 1 p^p 1 

As was mentioned in [l6| , for the self- force calculation for 
the electromagnetic case knowledge of the singular elec- 
tromagnetic potential is required, while for the self-force 
calculation for the gravitational case knowledge of the 
singular gravitational field is required. The correspond- 
ing equations are given in [l6l |. 

We argue here that the method for calculating the self- 
force that uses the singular and radiative fields is prefer- 
able to the one that uses the direct and tail fields. For 
the scalar, electromagnetic and gravitational cases, the 
singular field obeys the inhomogeneous Poisson, Maxwell 
and Einstein equations respectively and has support only 
outside the null cone of the particle. In that sense, knowl- 
edge of the entire past history of the particle is not re- 
quired, as is the case for the tail field, and the singular 
field needs to be calculated locally in the neighborhood of 
the point particle. Additionally, the radiative field obeys 
the corresponding homogeneous equation and thus is fi- 
nite and differentiable everywhere along the worldline of 
the particle. Consequently, this method for calculating 
the self-force does not present the difficulties of interpre- 
tation inherent in the method involving the direct and 
tail fields. 

A recent self-force calculation using the singular and 
radiative fields [IH showed that this method gives iden- 
tical results to the ones derived using the direct and tail 
fields. The results were also used [19( to predict the 
self-force effects on various orbits of scalar particles in a 
Schwarzschild background. Finally, a different approach 
to the self-force calculation described in [2(| uses the sin- 
gular field in order to extract the contribution of that 
field to the Weyl scalar Hfg (or ^|), subtract it from 
W T ct (Wj et ) resulting from the retarded field of the par- 
ticle and use the renormalized !^g en ( &T n ) to calculate 
the perturbations on the background spacetime. 

B. Outline 

In this paper we establish a method for calculating the 
singular field for monopole and dipole sources of scalar, 
electromagnetic and gravitational fields and we present 
a catalog of the expansions of the singular field for those 
sources. The coordinates that are used to calculate the 
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expansion of the singular field are the Thorne-Hartle- 
Zhang coordinates, (abbreviated as THZ coordinates). 
Those coordinates were initially introduced by Thornc 
and Hartle [2l[ and later extended by Zhang (22j. A 
short discussion of them is presented in Sec. [Til 

Expansions are given for the singular fields associated 
with all monopole and dipole sources for scalar and elec- 
tromagnetic fields and the monopole gravitational field 
in Sec. [TiTAl IhTbI and ITlTCl respectively. A short dis- 
cussion of the results is given in Sec. [TV] and future work 
is outlined in Sec. [VI 



II. THZ COORDINATES 

The Poisson, Maxwell and Einstein equations for the 
singular field assume a relatively simple form when writ- 
ten in a coordinate system in which the background 
spacetime looks as flat as possible. In the following, it is 
assumed that the particle is moving on a geodesic F in a 
vacuum background described by the metric g a b- Also, 
1Z is a representative length scale of the background ge- 
ometry, the smallest of the radius of curvature, the scale 
of inhomogeneities of the background and the time scale 
of curvature changes along the geodesic T. 

For a geodesic T, locally incrtial coordinates x a — 
(t, x, y, z) can be found so that, on the geodesic, the met- 
ric and its first derivatives coincide with the Minkowski 
metric and t measures the proper time along the geodesic 
[23l ]. Such a locally inertial coordinate system is not 
unique. The THZ coordinate system used here was first 
introduced by Thorne and Hartle [2l| and later extended 
by Zhang [22| and has meaning only locally, close to the 
worldline of the particle. Specifically it is assumed that 
the background metric close to the worldline of the par- 
ticle can be written as 



9ab — Vab + H a b 



(7) 



Vab + 2H ab + 3 H a b + °(^4~)> 



where r\ a b is the flat Minkowski metric in the THZ coor- 
dinates (t, x, y, z) and 



2 2,2,2 
p = x +y + z . 



(8) 



Also 



2 H ab dx a dx b = -£ ij x t x J (dt 2 + 8 k idx k dx l ) 
4 

+ 3 
_ 20 

~ 21 

5 

+ 21 



i^kpc[^^ x dtdx 



1 



dtdx K 



(9) 



■Ei^jpqB^ fcX^ X CpqilSj^X^ 1 



dx l dx 3 , 



3 H ab dx a dx b = - ~£ ijk x i x j x k (dt' 2 + 5 ni dx n dx l ) 
I ~^^-kpQ^^ i, 1 x dtdx 



I J - 

0(±-r)..dx l dxK 



(10) 



In this and the following, a, 6, c and d denote space- 
time indices. The indices i,j,k,l,n,m,p and q are spa- 
tial indices and, to the order up to which the calcula- 
tions are performed, they are raised and lowered by the 
3-dimensional flat space metric <5y . The dot denotes dif- 
ferentiation with respect to the time t along the geodesic. 
Also, Cijk is the 3-dimensional flat space antisymmetric 
Levi-Civita tensor. 

The tensors £ and B are spatial, symmetric and trace- 
free and their components are related to the Riemann 
tensor on the geodesic T by 



'R 



pqjt 
STF 



£ijk ~ [Vfc-Rtitj] 

Bijk = — [e.i pq ^ kR P qjt\ 



1 STF 



(11) 

(12) 

(13) 
(14) 



where STF means to take the symmetric, trace-free part 
with respect to the spatial indices i, j and k. The compo- 
nents £ij and Bij are of 0(1/1Z 2 ) and their time deriva- 
tives are of 0(1/1Z 3 ). The components £ij k and Bijk are 
also of 0(1/K 3 ). 

If H a b consists of the terms given in Eq. (|9|), the 
coordinates are second-order THZ coordinates and are 
well define d up to the addition of arbitrary functions of 
0(p 4 /TZ 3 ) 24]. Notice that in that case, only the Rie- 
mann tensor components and their time derivatives ap- 
pear in the metric. If H a b also includes the terms given 
in Eq. (fl"U|), the coordinates are third-order THZ coordi- 
nates and are well defined up to the addition of arbitrary 
functions of 0(p 5 /1Z 4 ) [24j. In this case, only spatial 
derivatives of the Riemann tensor components appear in 
the metric. 

The contributions to the metric components 2H a b that 
contain £jj and Bij are of 0(p 3 /lZ 3 ), as are the contri- 
butions to the metric components 3 H a b that contain £ij k 
and Bijk- Therefore, any calculation that includes the 
£ ij and Bij terms of 2 -Hah must also include the £ijk and 
Bijk terms of 3 H a b as well. 

Using the simple symmetry properties of the tensors £ 
and £>, the following relationships can be shown for their 
components and the spatial THZ coordinates (x, y, z): 



Ciik£ k i + £iki£ k j 



^ki£ k l = 

£ijkB 1 + eikiB k j - €jkiB k i — 



(15) 



and 



\^-ijk{^& n%l ^Zn*£ ) H~ €ikl& 'n%j £jkl& n>£i\& % 

[e ljk {B k n xi - Bi n x k ) + e lk iB k n x J - e jH B k n Xi]x n x l = 0. 

(16) 
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These relationships are used in Sec. IIIII where the singu- 
lar fields for different sources are calculated, in order to 
simplify the expressions for those fields. 

Before we proceed to the calculation of the singular 
fields, we need to make an important point on the no- 
tation and language that we use. In the following, the 
subscripts (or superscripts) (0), (2), (3) and (4) are used 
to indicate the order of significance of each term or com- 
ponent, for the singular fields, the covariant derivatives 
with respect to the THZ metric and the components of 
the Riemann tensor that show up in the Poisson, Maxwell 
and Einstein equations below. The subscript (0) refers to 
the most dominant contribution (resulting from the part 
of the THZ metric that is of order 1). The subscripts (2) 
and (3) refer to the next two more significant corrections 
(resulting from the parts of the THZ metric that are of 
order (p 2 /lZ 2 ) and (p 3 /lZ 3 ) respectively), which are cal- 
culated for the singular fields. The subscript (4) refers to 
the next correction (resulting from the (p 4 /1Z 4 ) part of 
the metric). The subscript (1) is not used due to the fact 
that there are no (p/TZ) terms in the metric. The terms 
"first correction" , "second correction" and "third correc- 
tion" refer to the correction coming from the 0(p 2 /1Z 2 ), 
0{p 3 /1Z 3 ) and 0(p 4 /lZ 4 ) parts of the metric respectively. 
The fact that they are labeled with indices (2), (3) and 
(4) (instead of the intuitive choice (1), (2) and (3), respec- 
tively) should not cause any confusion, since the reason 
for it is justified. 



III. CALCULATION OF THE SINGULAR 
FIELDS 

The singular fields for a variety of sources are calcu- 
lated below. The calculations were done using the pro- 
gram GrTensor [HI running under Maple. 

It should be noted that in all cases it is assumed that 
the characteristic of the moving source that is responsi- 
ble for its own field (i.e. charge, dipole moment, mass) 
is small compared to the background curvature. That 
assumption allows us to safely consider the effects of the 
self-force as a small perturbation on the motion of the 
particle. 

The derivation of each singular field follows a similar 
pattern. To illustrate that pattern we use the simple ex- 
ample of the scalar monopole field. In THZ coordinates 
the monopole scalar field is assumed to have q/p as the 
leading term in the expansion in powers of the distance 
away from the point source. This singular behavior cor- 
rectly accounts for the 5-function source at the origin. 
However, the V 2 operator in curved spacetime, in THZ 
coordinates, differs from its flat space counterpart by 
terms involving the £ij and £>y multipole moments; these 
have dimensions of l/(length) 2 and are of 0{l/1Z 2 ). It is 
obvious that the first correction to the singular monopole 
scalar field must behave as g(£y or Bij)x l x 3 / p in order to 
have the proper index structure and dimension. The cor- 
rection must have no free indices, and the two indices on 



either of the the tracefree £ij or £>y require precisely the 
x % and x 3 . The power of p in the denominator is deter- 
mined by the need to have the correct overall dimension. 
With a similar argument it follows that the succeeding 
correction to the singular monopole field should generally 
be of the form q(£ijk or Bijk)x' l x 3 x k / p. 

This picture of the expansion is consistent with the 
general description of the Hadamard expansion [26| of 
the Green's function in curved spacetime for the singu- 
lar field 0, [l6| . The issue that remains is to verify that 
this expansion does not surreptitiously include some of 
what ought to be considered the "regular remainder" R- 
part of the retarded filed, which is responsible for the 
self-force. Assume, for the moment, that the expansion 
of the correct singular monopole field does include a reg- 
ular piece; such a regular piece must be a homogeneous 
solution of the field equations in the vicinity of the par- 
ticle because the (5-function source is already accounted 
for in the expansion which was just described. Thus, the 
general form of such a correction must be expandable in 
THZ coordinates about the particle, and must behave as 

constant + linear term + quadratic + cubic + . . . . (17) 

The first derivative of the constant term vanishes at the 
origin, and therefore could exert no force on the parti- 
cle. Similarly the first derivatives of the quadratic, cubic 
and higher order terms also vanish at the location of the 
particle and cannot contribute to the self-force. Only the 
first derivative of the linear term remains. However, the 
linear term ought to be describable only in terms of the 
geometry and should be expressible as 

q(£ij or Bij) x (something linear in x l ). (18) 

But there is no "something linear" in x % which has the 
proper dimensions and index structure. Thus, while the 
singular monopole field constructed above might, in prin- 
ciple, include an extra "homogeneous part," the extra 
part could not be linear in x % and could not affect the ac- 
tual self-force on the particle. A similar argument can be 
used to rule out the presence of at least regular quadratic 
and cubic terms as well. 

It will be seen in the following that the order up to 
which we have to calculate the singular fields (in order to 
be able to calculate the self-force) depends on the particle 
in question. In some cases, it is sufficient to calculate 
the first correction, namely that which results from the 
0(p 2 /1Z 2 ) parts of the THZ metric. In the other cases, 
we must also know the second correction that comes from 
the 0(p 3 /lZ 3 ) parts of the metric. For completeness, we 
calculate all these corrections for all the sources that we 
examine and we explicitly state in which cases that is not 
strictly necessary. 

In the cases in which calculating the second correction 
is not strictly necessary, we can still significantly bene- 
fit from such a calculation. Firstly, the more accurate 
the approximation of the singular field, the more differ- 
entiable the radiative field. For example, as explained in 
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[18| , if the approximation to the scalar singular held is in 
error by a C n function, the approximation to the scalar 
radiative held is no more differentiable than C n and the 
approximation to the self-force (obtained by differenti- 
ating the radiative held) is no more differentiable than 
C n_1 . Additionally, according to the standard method 
for calculating the self- force [7|,one needs to do a spher- 
ical harmonic decomposition for the retarded and singu- 
lar fields, subtract the 1-mode contribution of each one 
and finally sum over all 1-modes. As was seen in [l8[ 
and in [27| , the convergence of the sum over the spheri- 
cal harmonics can be sped up by including the spherical 
harmonic decomposition of higher order corrections. 



involves £y/- and Bijk as well as the time-derivatives of 
£ij and Bij. From this equation it can be inferred that 



(3) 



[O(*(o)) x 0(£a)]. 



(25) 



Notice that \I> S 2) 



because it shows up in the term V/ 2 ) * (2) • 
[O(*f 0) ) 



(3)' 

a term of order 



is not included in the equation for tyf, 
2 mS 

(2) 

) x 0(p 4 /1Z 4 )], which must be included in the 
calculation of the third correction. 

The source for the third correction contains both \&^L 
and ^( 2 )" Specifically, the differential equation derived 
from Eq. (|20[) by keeping all metric terms of 0(p 4 /7Z 4 ) is 



Scalar Field 



V (0)*(4) 



72 ,t/S 1 w2 <qjS 



(4)^(0) 







(26) 



Assume that a particle moving on a background 
geodesic creates a scalar source g. The scalar field gen- 
erated by the particle can be expanded as 



\J7 



(2) *(3) 

and it obeys Poisson's equation 



(19) 



V 2 * 3 = V 2 0+2+3+ . 



s 

(0) 



(2) 



(3) 



.) = -iirg. 
(20) 

Regardless of the form of g, the leading term in the ex- 
pansion is the scalar field generated by the particle when 
it is stationary at the origin of a Cartesian coordinate 
system and it obeys the lowest order differential equa- 
tion derived from Eq. (|2"0|) 



V (0)*(0) 



-AltQ. 



(21) 



The first correction obeys the differential equation de- 
rived from Eq. (|20|) by keeping only the lowest-order 
terms, namely 



Vfo)*(2)+V 2 2) *f ) 







(22) 



which means that the V 2 2 -j acting on the zeroth-order 
part of the field is equal to minus the source term in 
the scalar differential equation for the first correction to 
the field. In general, the source term in this equation 
is expected to contain the components of the tensors E 
and IB and to give the first correction coming from the 
particle's motion on the background geodesic. From this 
equation it can be inferred that 



^ 2) ~[O(* s 0) )x 0(^)1 



(23) 



The second correction obeys the differential equation 
derived from Eq. (|20p by keeping the 0(p 3 /1Z 3 ) terms of 
the metric, namely 



V 2 o)*(3)+V 2 3 )*fo) 



0. 



(24) 



The source term in this case is identically equal to minus 
the 0(p 3 /7tL 3 )-part of V 2 acting on and in principle 



in which the second and third terms are effectively the 
source terms. That gives that the third correction is 



(4) 



[°(*(0)) X °(^4)]- 



(27) 



1. Monopole Field 

Assume that the particle in question is carrying a 
scalar charge q and is moving on a geodesic. The scalar 
field generated by that particle is the scalar monopole 
field. The detailed calculation of that held in the THZ 
coordinates was presented in [l8j and the result is given 
here for completeness. 

The zeroth-order term is the Coulomb-like potential 
(<?//?)■ Both the hrst and second order corrections can be 
shown to be equal to due to symmetries, so the result 
is 



(28) 



Disregarding the fact that the first and second correc- 
tions to the held are equal to zero, Eq. (f23|) and (f25|) im- 
ply that one would only need to calculate up to the first 
correction of the scalar monopole singular field in order 
to be able to calculate the self-force on the scalar charge. 
That is because the second correction, were it not zero, 
would be of 0(p 2 /1Z 3 ) which would give no contribution 
to the self-force on the particle after the first derivative 
followed by the limit to the location of the particle were 
taken. 



2. Dipole Field 

Assume that the particle is carrying a scalar dipole 
moment and is moving on a geodesic of the given back- 
ground. The dipole moment is assumed to have a ran- 
dom orientation and its THZ components are denoted as 
K a = (Q,K x ,K y ,K z ). 
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The leading term in the scalar field series expansion is 
the scalar field generated by a dipole that is stationary 
at the origin of a Cartesian coordinate system and is 



(0) 



Kjx 1 
P 3 



(29) 



For the first correction, explicit evaluation of Eq. 



shows that V 2 2 )^( ) = ^> wmcn results in 



' dx 2 



dy 2 



d 2 



o. 



(30) 



so the first correction to the field can be set equal to a 
constant, or, for self-force calculations, equal to 0. 

Eq. gives that the second correction to the singular 
field obeys a differential equation that relates the second 
derivatives of to terms of the form 



KB 



where the dots denote appropriately contracted indices 
for each term. The tensors £ijk or Bijk do not show up 
in the equation. The fact that there is a large number of 
terms in the differential equation makes it very inconve- 
nient to present the equation explicitly here. In order to 
give the reader an idea of what the equation looks like, 
we explicitly write a few terms of it [28j | , namely 



& 2 



d 2 



d 2 



dx 2 dy 2 dz 2> ^ 



K X B 



x z y z z z 



= 0. (31) 



Using arguments similar to the ones mentioned at the 
beginning of Sec. (jlllj) for the scalar monopole field we 



conclude that must be proportional to 



(3) 



Ae jkl KW k . 



I i 
XX 



where A is a constant. Substituting this expression into 
the differential equation for gives that A = —1/7, or 
that 



(32) 



Eq. |27J gives that the third correction is of 0{p 2 /K 4 ). 
This correction is not necessary for self-force calculations. 
Differentiation of it would give terms proportional to the 
coordinates, which would give zero after the limit to the 
location of the particle is taken. 

Finally, the singular scalar field of a dipole moving on 
a geodesic is equal to 



K,x* 



1 K^B k iX l x l 
ejkl- 



0( 



K 4 ' 



(33) 



P ' P 

In this case, knowledge of the second correction is nec- 
essary in order to calculate the self-force. Differentiating 
the second correction gives a term of order O(p /7Z 3 ) 
which is non-zero when the limit to the location of the 
particle is taken. 



B. Electromagnetic Potential 

Assume that a particle carries an electric charge, an 
electric dipole moment or a magnetic dipole moment. 
The source associated with it is denoted by J a . In this 
case, the self- force can be calculated using the singular 
electromagnetic potential, which can be expanded as 



4° 



^S(2) 



A a 



(34) 



In the vacuum background the Ricci tensor is equal to 
zero. That means that the electromagnetic potential 
obeys the Maxwell's equations 



X7 2 A a 

v (0+2+3+.. .)^ 1 S(0+2+3 



+3+...) 



-47rJ a . 



(35) 



Regardless of the exact form of the source, the leading 
term in the expansion (|34[) is the electromagnetic po- 
tential that would be generated by the particle if it were 
stationary at the origin of a Cartesian coordinate system. 
It obeys the equation 



V (0)^3(0) 



-47rJ Q . 



(36) 



The first correction obeys the differential equation de- 
rived from Eq. (|35|) . namely 



V 2 A a 
V (0) /1 S(2) 



^ (2)^13(0) 



0. 



(37) 



The 0(p 2 /lZ 2 )-pa,rt of V 2 acting on the zeroth-order elec- 
tromagnetic potential is the source term for the first cor- 
rection. From this equation it can also be inferred that 
the order of the first correction to the electromagnetic 
potential is 



(38) 



According to Eq. (|35|) . the second correction obeys the 
differential equation 



V7 2 A a 

V (0)^S(3) 



^(3)^S(0) 







(39) 



meaning that the 0(/? 3 /7?. 3 )-part of V 2 acting on the ze- 
roth order electromagnetic potential acts as the source 
term in the differential equation for A^ 3 y As in the 
scalar case, ^4g( 2 ) does not show up in this equation, but 
will show up in the equation for the third correction. The 
equation for A^^ indicates that 



A a 



(40) 



The equation for the third correction as derived from 
Eq. (0SD is 



V 2 A a 

V (0)^S(4) 



^(2)^S(2) 



^(4)^|(0) 







which gives the order of the third correction 



A a 

^S(4) 



[0(A a m )xO(^)}. 



(41) 



(42) 
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1. Monopole Potential 

If the particle is endowed with an electric charge q, 
Eq. (f3"6")l gives that the zeroth order potential is the 
Coulomb electromagnetic potential 



A a m = [1 0, 0, 
P 



Substituting the THZ components of A% 



A a 

A S{2) 



S(2) 



(^■S(2)' yl S(2))^ l S(2 



(43) 



(44) 



into Eq. (|37|) results in four differential equations, one 
for each one of these components. Each equation relates 
the sum of second derivatives of a component to a sum 
of terms of the form: 



and substituting into Eq. ([39]) gives a differential equa- 
tion for each component of ^g( 3 )- Each equation relates 
a sum of second derivatives of each component to terms 
of the form 



./ « !• ■ I « < - « * ■ 



for the f-component, and 

q£ ^ — and qB ^ — 

p J p 6 

for the spatial components. The t-component in this case 
is expected to be proportional to 



q£. 



x l x 3 x k 



ijk' 



q£. 



X^X^X^X^ 



for the i-component and 



qB.. 



for the spatial components, where the dots again denote 
appropriate indices. 

The form of the solution can be predicted based on the 
dimensionality and the index structure as before. One 
must take into account the fact that the ^-component 
must have no free indices and that each spatial compo- 
nent must have one free index. Thus it can be derived 
that the solution should be proportional to 



q£ t 



for the f-component and to 



qe p ijB\ 



x k x 3 
P 



for the p-spatial component, each term multiplied by 
an appropriate constant. Substituting these expressions 
into the four differential equations gives simple algebraic 
equations for these constants, which can be easily solved 
to give the components of the first correction 



4* 

A S{2) 



A p 

^S(2) 



XX 



x 3 x 



\qz p nB\ 

z p 



(45) 



The second correction to the electromagnetic potential 
comes from the part of the V 2 that is of O(^j) acting 
on the zeroth-order electromagnetic potential. Assuming 
that the THZ components of A§^ are 



^S(3) — (As(3)> ^S(3)>^S(3V ^1(3)) 



(46) 



and the p-spatial component is expected to be a sum of 
the terms 



q£ P iX l p, q£ l 



qe p tj B lj k x k p, qePijB'kl- 



x 3 ' x k x l 



qtijkB lk i 



X p X 3 X l 

p 



each term multiplied by an appropriate constant. Sub- 
stituting these expressions into the four differential equa- 
tions results in algebraic equations for the constants. Fi- 
nally, the second correction to the electromagnetic po- 
tential is 



^s(3) — ~ ~Zq£i]k 



x l x 3 x k 



A p 

^S(3) 



1 ~ n j 7 • x^x x 
-q£P lX p + - q£lJ —- 

-qe i:j B u . 

9 p 



(47) 



According to Eq. (|42|) the order of the third correction 
to the electromagnetic potential of a charge is 



Al (4) ~ 0( 



ft 4 ' 



(48) 



and it is not necessary for self-force calculations. 

Finally, the components of the singular electromag- 
netic potential of a charge q that is moving on a geodesic 
are equal to 



M -1 Lf x%xJ 1 

p z 



'.I 



^q£ijk 



x z x 3 x k 



p 3 

0(|j) 



x 3 x k 



—q£ P iX l p + —q£ i: j 
lo 



1 

Ys 

l&vBL— 



x p x l x 3 



(49) 



x 3 x k x l 



0( 



K 4 ' 
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In this case, the second correction is not necessary for 
the calculation of the self-force. After differentiation, 
that correction will give terms proportional to the co- 
ordinates which, after taking the limit to the location of 
the particle, will give zero. 



algebraic equations for the multiplicative constants. The 
final expressions for the components of ^4g( 2 ) are 



A 



S(2) 



1 x % x 3 x k 

nH^Ok -3 , 



2. Electric Dipole Potential 

Assume that a particle carrying an electric dipole mo- 
ment is moving on a geodesic T. The dipole moment is 
assumed to point at some random direction and its THZ 
components are q a = (0, q x ,q y ,q z ). 

According to Eq. (|36[) the leading term in the expan- 
sion of the singular electromagnetic potential is the elec- 
tromagnetic potential generated by an electric dipole that 
is stationary at the origin of the Cartesian coordinate sys- 
tem 



1 „ „„■ x°x k x l 



The second correction to the singular electromagnetic 
potential obeys Eq. ([55]) . By substituting its THZ com- 
ponents 

^S(3) = (^S(3)'^S(3)'^S(3)'^S(3)) (53) 



( 'yr- »• 0,0 



The THZ components of the first correction are 



into (|39|) we get a set of four second-order differential 

(50) 

equations for each one of those components. Each equa- 
tion relates the second derivatives of one component to 
terms of the form 



^S(2) — (^S(2)'^S(2)>^S(2)'^S(2))- (51) 

If substituted into Eq. (f3"T)) . the differential equation for 
^S(2) ^ s broken into a set of four second-order differential 
equations for these components. Each differential equa- 
tion relates the sum of second derivatives of a component 
to a sum of terms of the form 



tje ' X> ' *X ' X " tlj ' 



for the f-component, 



q.B. 



rp ' rp ' rp ■ rp ■ rp ■ 



rp • rp ' rp ' rp • rp • rp • rp ' rp • rp • rp • rp • rp • 

q.B.. 7 ,q.£... 7 

P 7 P 7 



for the ^-component, and to terms of the form 



rp ' rp • rp ' rp ' rp ' rp ' rp ■ rp ' 

q.£.. — s — ,q.B.. 



p° 



for each spatial component. 

In this case the solution must be a sum of the terms 



e i3k q l B 3 i- 



k I 

XX 



for the spatial components. 

In this case we can deduce that the solution should be 
equal to a sum of the terms 

rpi rpj rpk ,y,J 

W Us Uj % O 

qitjk — ; — , q tij — 
p 6 P 

for the i-component and a sum of the terms 



x 
P 

x 3 X X 1 



x 3 x x l 



e^rfB 3 ^, e^rfBu^j^, e p i jq k B i k — , 



e p lJ q k B l l ^-f-, W^y. eaktf&i , 



X P X X 



for the p-spatial component, each term multiplied by an 
appropriate constant so that the differential equations 
are satisfied. Additionally, using Eq. (fT5|) and (fH)|) , the 
last two terms that are expected to show up in the so- 
lution for the p-component can be eliminated in favor of 
the remaining four. Substituting the sums of the remain- 
ing terms into the four differential equations gives simple 



rpi rpj rpk rp^ rpj rp^ 



qi£jkl 3 ; q l £ijk 

P A P 



for the f-component, and a sum of terms 



X^ X^ ' X ° X^ 

q p £ij , q j £ p l , q l £ P ip 

P P 

rpP rpj rpP rp . rpj rp^ 

q £ij——, q £jk — j 3 — 



k I 

XX 



pj rpk ™,/ ,-pTTL 



e p ij q l B } k i—i e p lj q l B Mm - 



P 



-, e p ij q k B i M 



r j n,,l 



p kni X 3 X k X X 

t ijq B i m -jg , e ijk q B J i m _ . 



x p x k x l x m 



for the p-spatial component, each term multiplied by an 
appropriate constant so that the differential equations 
are satisfied. We substitute these expressions into the 
differential equations and solve the resulting systems of 
algebraic equations for those constants. Thus, the second 
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correction to the electromagnetic potential is 



1 



As (3) = - 



■a X k X l 1 



x i x^x x 



p 6 p J 



S (3) " lg 



p 18' 



1 , • x p x j 7 x p XiX j x k 

—q £ij h — q t jk 5 

18 p 18 p A 



+ -e p l3 q k B\ rn 
9 p 



x j x k x l x m 2 
3 + -e l3k q l B n r 



According to Eq. (|4"2")l the third correction is 

(4) - O(^) 



(54) 



(55) 



which is clearly not necessary for self- force calculations, 
because taking the first derivative and the limit to the 
location of the particle gives zero. 

Finally, the singular electromagnetic potential for an 
electric dipole moving on a geodesic is equal to 

. t x % 1 x % x^x k 

J ,-~„- x k x l 1 „ x l x^x k x l ^, p 



7 



€ijkq t B : 'i— — -^qi£ju j 



P 6 



1 



x l 1 



s 2 JH p 2 



/0 



x^x k x l 



1 



18^ 



+ -e p l3 q k B\ m - 



18^ 



h 1 

XX 



(56) 



In this case the second correction is necessary for the 
self-force calculation. That correction is proportional to 
the first power of the coordinates, differentiation of which 
gives a constant term that does not vanish when the limit 
to the location of the particle is taken. 



3. Magnetic Dipole Potential 

Assume that a particle with a given magnetization is 
moving on the geodesic T. The magnetization m a is as- 
sumed to point at some random direction and its THZ 
components are m a = (0, m x , m y , m z ). 

Eq. (|36|) gives that the leading term in the expansion of 
the electromagnetic potential is the potential generated 



by a magnetic dipole that is stationary at the origin of a 
Cartesian coordinate system. Its THZ components are 



S(0) — V"' c y""* n3 > e ij m ^ 3 j c y" ^ 3 



! 'V'''^- (57) 



The first correction ^ has THZ components: 

^S(2) = (^S(2))^S(2)>^S(2)'^S(2))- (58) 

When those componenets are substituted into Eq. (pT7|) 
the result is a set of four second-order differential equa- 
tions, one for each one of those four components. Each 
equation relates a sum of the second derivatives of a com- 
ponent to a source term that consists of terms of the form 

„ x x x- 
m B — = — 



P° 

for the f-component, and terms of the form 



for the spatial components, where the dots denote appro- 
priate indices. 

The form of the solution can be predicted as previously. 
In this case, the i-component is expected to be a sum of 
the terms 



m,iB l j 



x^ m i BjkX t x 3 x k 

> 3 
P p A 



and the p-spatial component a sum of the terms 



x^x k x l 



e p ij m l £ : > k — , e p ijm l £ kl T 

P P 3 

,>■,,„<>■£■ i-" 1 -'' 1 - 1 '' 



x k x l x p 



, , tijkm % Z pl — , e ijk m l £ :l i ~ 

p 6 p p A 



with appropriate constants multiplying each term so that 
the differential equations are satisfied. Using Eq. (|T5|) 
and (fT6|) . the first and last terms expected to show up 
in the sum for the p-component can be eliminated, since 
they can be expressed as linear combinations of the re- 
maining four terms. Substituting these expressions into 
the differential equations gives simple systems of alge- 
braic equations for those constants. The result is that 
the first correction to the electromagnetic potential has 
components: 



.4 



1 



\rriiBjk 



x l x j x k 2 



S(2) -g"'»^Jft n3 



m l Bi 



^S(2) — £ P ij ml £kl 



X^X k X l 1 



,m k Z'' 



x^x k x l 



2 1 n 3 



(59) 



1 . . T fe 

- -e iik m > £ p: >—. 

2 3 p 

The second correction to the magnetic dipole potential 
can be calculated using Eq. ([39)) . Substituting the THZ 
components of that correction: 



^S(3) — (^S(3)'^S(3)'^S(3)'^S(3)) 



(60) 
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to Eq. ((39)) we get four differential equations, one for 
each component, which relate second derivatives of those 
components to terms of the form 

* X X lC tjC Xj X X tC 

m.t.. z ,m.o... = 



p° p° 
for the i-component and terms of the form 

X X X- X X X _ tjj *jC ^jC X- X> *jC 

mo = ,mb = 



p> p< 

for the spatial components. The terms expected to show 
up in the solution can be constructed as before. We ex- 
pect a sum of the terms 



e i j k m' l £ 3 i 



k I 

X K X l 



x 3 x k 



m l B ljk , miBi 



x l x 3 x k x l 



P P 3 

for the ^-component and a sum of the terms 



xx- 



X Xj 



wiPlBij—^— , m l B p lP , m J B p l ^- 

■ ■ x p x 3 1 

rriiB 1 j , -^miBnkX p x l x :) x k 

P P 3 

X X X^ X 

e p in m l £ 3 k i , e p ij m k £ l M - 



k I 

eijkm l £ pj i , e i jkm t £ j in 



x p x k x l x n 



C P ijTTt £kln 



x 3 x k x l x n 



P' 



; € P ijm n £ kl 



x 3 x k x l x n 



for the p-spatial component, each term multiplied by an 
appropriate constant. Substituting those sums into the 
differential equations, we get systems of algebraic equa- 
tions for those constants. Solving those we get that the 
second correction to the electromagnetic potential is 



.4 



1 



S(3) — ytijklri&l 
1 



k 1 
XX 



- -m % Bijk 1- -mioiik- 

3 P 9 p 3 

37 . 95 . 5 . rPr-? 

= ^Lm p B in — + ^m^ip + i-rrfBij — 



x 3 x k 1 „ x l x 3 x k x l 



-miBi 



s (3) 252 



J '.i 



p 126 



63 

x p x l x 3 x k 



— m B p j -^rruBjk 5 

63 p 36 p 6 

-e p ijm l £kin 5— e p li m k £ l ki — 

6 J p 3 6 3 p 



X 3 x l 



1 



x p x k x l x n 



-e ljk m l £ 3 i n o 

6 p A 



(61) 



The order of the third correction can be predicted as 
usual, by Eq. ([IT]) and it is 



4(4) - 0(£l)- 



(62) 



These correction terms are not needed for self-force calcu- 
lations. As in the case of the electric dipole electromag- 
netic potential, differentiating these terms would result 
in terms that are proportional to the cooridinates and 
taking the limit to the location of the particle gives zero. 

Finally, the singular electromagnetic potential gener- 
ated by a magnetic dipole moving on a geodesic Y has 
components equal to 



I rj.i rpj ryk O 

At W 1 VJ 

A s = -rrnBik 5 £> y — 

3 p 



i s - -nliDjk— -3- 



1 i6i x k x l 
+ -j^ijkm £ 3 i—— 

1 , x 3 x k 1 

- -m Bijk h -miB l jk 

3 P 9 



x % x 3 x k x l 



p _ „ mV 



(71, 



x^ X X 



m,k£ % lX 3 X k X l 1 



2 y p 



p 126 



63 3 
1 



63 J 



1 



x p x l x 3 x k 



36"" ifc p 3 



e p il m l £', 



6' " 
1 
6 



- -^i]k1Tl l £ 3 h 



x 3 x k x l x n 1 
p 3 6 
x p x k x l x n 



t p ^m k £ l kl 



X J X' 



(63) 



It is clear that in the case of the magnetic dipole , it is 
necessary to know the second correction to the singular 
electromagnetic potential in order to be able to calculate 
the self-force. 



C. Gravitational Field 

We proceed by looking into the more interesting (and 
computationally more complicated) case of a particle 
moving on a background geodesic Y and producing a 
gravitational field. The singular gravitational field can 
be expanded as 



n Sab — n s ab ' n S ab ' n S ab 



(64) 



For the calculation of the gravitational singular field it is 
more convenient to use the trace- reversed version of h^at 
which is defined as 



We work in the harmonic gauge 
V a h Sab = 0. 



(65) 



(66) 
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The singular field obeys the linearized Einstein equations 
which, in the harmonic gauge, become 



V 



(0+2+3+.. .)"S 



"(0+2+3+...) 



2R 



(0+2+3+.. .)a b rig 



ab 

c dT(0+2+3+...) 



cd 



-l&TlTab- 



(67) 



In Eq. ([67]) T a b is the stress-energy tensor associated with 
the particle. 

Regardless of the exact form of the stress-energy ten- 
sor, the leading term in the expansion of the gravitational 
field is the field that the particle would produce if it were 
stationary at the origin of a Cartesian coordinate system 
and it obeys the equation 

V( )4°L + 2i? (0)a c 6 rf 4 0) «2 = -16^T afc . (68) 

The first correction obeys the differential equation de- 
rived from (|67p by keeping only terms that come from 
the 0(p 2 /7Z 2 ) part of the metric, namely 



T7 2 I o p c dT(2) _ 

V (0) ft S ab ' ZU (0)a b n s c d — 

V7 2 OB c dT(0) 

— V (2) n 'S ab ~ ZJrt -(2)a b fig cd- 



(69) 



The V 2 and Riemann tensors of 0(p 2 /1Z 2 ) acting on the 
zcroth order gravitational field act as the source terms in 
the differential equation of the first correction. Eq. ([69]) 
also implies that the first correction is 



h^ab- [0(4°' afc ) X 0{^)} 



(0) 



n 2 ' 



(70) 



Similarly, the second correction obeys the differential 
equation derived from Eq. (|67[) by keeping only terms 
that come from the 0(p 3 /1Z 3 ) part of the metric, namely 



V (0) ft S afc + ^- K (0)a b n s c d — 

V7 2 OB c dr(0) 

— V (3)i*g ab — z -n-(3)a b cd- 

Eq. (171]) also implies that the second correction is 



(71) 



4 3) ^[OC»MxO(^)] 



(o) 



(72) 



The order of the third correction to the gravitational 
field can also be predicted, based on the differential equa- 
tion that could potentially yield that correction. Specif- 
ically, Eq. (frJT)) gives 



V7 2 i 9B c d h^ — 

V (0) ft S ab + Zn (0)a b fig cd — 

vy2 T(0) 9R c dt(0) 

— V^/lg a b — /-ft(4)a b "s cd 

V/2 7(2) 9R c dr(2) 

— V(2)"'S ab ~ z/t (2)a b IT'S cd 

which implies that the third correction is of order 



(73) 



h { s 4) ab~[0(h^ a b)xO(^)] 



(0) 



ft 4 " 



(74) 



1. Massive Particle 

Assume that the particle that generates the gravita- 
tional field is a massive point particle of mass m. The 
calculation for the first correction to the gravitational 
singular field for this case was first shown in [24[ , where 
a different coordinate system was used. Here we perform 
the calculation of the singular field in the THZ coordinate 
system, including up to the second correction to it. 

The leading term in the expansion of the singular grav- 
itational field is the gravitational field generated by a 
massive particle that is stationary at the origin of the 
Cartesian coordinate system. We consider only the terms 
that are of first-order in the mass m. Thus 



h {0) , 

"S a-b 



rn 



/l 000 
10 
10 

Vo o o i 



(75) 



The first correction to this gravitational field obeys 
Eq. (|69)l and must be a symmetric tensor. Substituting 
its components into Eq. (|6"9"]) results in ten second-order 
differential equations. There is one independent set of 
four differential equations for the four diagonal compo- 
nents, each equation containing all four diagonal com- 
ponents. There is also one differential equation for each 
of the t — i components and one differential equation for 
each of the spatial i — j components with i ^ j. In each 
equation a sum of second derivatives of components is 
related to a sum of terms of the form 



for the t — t and p — q components and a sum of terms 



mB 



for the t — p components. 

Solving the differential equations in this case is more 
complicated than in the previous cases, mainly because of 
the fact that four of them involve all diagonal components 
rather than only one of them. However, the process that 
was described at the beginning of Sec. IIIII for predicting 
the form of the solution is applicable in this case as well. 
It is clear that the t — t component must have no free 
indices, each of the t — p components must have one free 
index and each purely spatial component must have two 
free indices. We conclude that the t — t component of the 
first correction to the singular field must be proportional 
to 



X X J 



P 



each t — p component must be proportional to 

■ x^x k 



me vii B l k- 



P 



12 



and each purely spatial p - 
of the terms 



q component must be a sum 



and the p — q spatial components are expected to be a 
sum of the terms 



mSpqp, m£p 



rr> T rv* nr>J 

mtij \ 



TflTjpq&i 



p° p 

with appropriate constants multiplying each term so that 
the differential equations are satisfied. Eq. (fT5| and (fT6|) 
can be used to eliminate the second and third terms in 
favor of the remaining two, for the expression for the 
p — q spatial components. Substituting into the differen- 
tial equations we get simple algebraic equations for the 
multiplicative constants. The result is that the first cor- 
rection to the singular field has components 



S (2) 

^ s tt 



,(2) 



2m£i 



An-! 



me plj B l k 



x^x 



(76) 



x^x 



l J /'•' 



XX- 



TjpqUlE.ijk ^ , TflSpgiX /), UlS^j^pXq^ 



each term multiplied by an appropriate factor so that 
the differential equations are satisfied. The paranthcscs 
around two indices denote the usual symmatrization con- 
vention with respect to those indices, for example 



j(p^q)-^ ^ ci \^ijp^q ~f~ £-ijqXp)x X^ . 



■,(2) 
'S PI 



-Am£p q p - 2mrj P q£i 



XX- 



The second correction is also a symmetric tensor 
and must obey the differential equations derived from 
Eq. (fTTj) . Substituting the components of the second cor- 
rection results in ten differential equations. As in the case 
of the first correction, four of these equations involve sec- 
ond derivatives of the four diagonal components. Each 
one of the other six differential equations contains the 
second derivatives of the six non-diagonal components of 
the second correction. The four equations for the diag- 
onal components and the three equations for the non- 
diagonal purely spatial components relate derivatives of 
those components to terms of the form 



iB 



m£ 



P P 

The equations for the t— p components relate their deriva- 
tives to terms of the form 



m£ 



raB 



P P 
We can predict the form of the solution as previously. 
The t — t component must have no free indices. It can be 
easily deduced that all possible terms that contain B and 
satisfy that requirement are identically equal to zero, so 
the t — t component is expected to contain only the £ijk 
tensor and be proportional to 



m£- 



x l x^x h 



ijk' 



P 



The t 
terms 



p components are expected to be a sum of the 



m£pix l p, m£ij 



ry*1 ryj qr> 
Jb Jb Jb Tj 



me vii B 



x J x X 



pij & hi ' 



Substituting these expressions into the differential equa- 
tions gives simple algebraic equations for the multiplica- 
tive constants, which can be easily solved. The second 
correction is then equal to 



,(3) 



,0) 

*S tp 



l(3) 
'S PI 



- g TFt&ij k 

10 



x l x^x k 



38 



63 p ' 63 J , 



X J X k X 



(77) 



— —Tnc r) jil3''fri 

9 P3 p 

10 • , , 10 - % xix k 

■—me lj{p B q) x J p- —me lj{p x q) B k—^~ 

2 x l x°x k i 

T^lpqm&ijk ZUlOpq-lX p. 



The order of the next correction can be predicted by 
Eq. (J73D and it is 



*S(4) 



(78) 



Evaluating these terms is clearly not necessary for self- 
force calculations because they will result in zero con- 
tribution after the first derivative and the limit to the 
location of the particle is taken. 

Finally, the singular gravitational field of a point par- 
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ticle of mass m moving on a geodesic is equal to 

C^1>J 2 f>J 'V*'* /l3 



point r 



2m 
P 



x l x^x 



p 

2 - x^x k 

hstp = - -^me pl jB t k 

6 p 

10 • „■ 38 



P 



// S // =— + liaEij^— + |m% fe "" 4- 



lu A j r> } A 

— mt„jX p men — 

63 p H 63 3 



-me pij B l k r 
2m 



X^X X 



hspq =r lpq — 4m£ pq p - 2mr] pq £ ij —— 



x 3 x h 



10 R < 7 10 Ri 

—me i3[p B q) x J p- —me ijip x q) B fc — 



x l x^x 



2m£ pqi x l p + 0(^). 



(79) 



Direct substitution into Eq. (|6"6"|) shows that this expres- 
sion for the singular field is consistent with the harmonic 
gauge condition. 

It is clear that, in this case, the second correction to 
the singular gravitational field is not necessary for the 
calculation of the self-force on the massive particle. Af- 
ter differentiation, that correction will give terms propor- 
tional to the coordinates which, after taking the limit to 
the location of the particle, will give zero. 



IV. DISCUSSION 

There is a very obvious categorization of the particles 
mentioned above, based on the order up to which it is nec- 
essary to calculate the singular fields in order to be able 
to calculate the self-force. For the particles that carry no 
intrinsic spin, namely the scalar monopole particle, the 
charged particle and the non-spinning massive particle, it 
is sufficient to calculate the singular fields/potentials up 
to, and including, first corrections. For the particles that 
do carry intrinsic spin, namely the scalar dipole particle 
and the particle carrying an electric or a magnetic dipole 
moment, it is necessary to include the second correction 
to the singular fields, in order to be able to calculate the 
self-force. 

That can be explained quite easily by considering the 
fields generated by each particle and we use here the sim- 
ple example of the scalar fields. Let us compare the scalar 
monopole field at a point r 



ijj(f) = 



<1 



\r-r \ 



(80) 



(81) 



of a scalar charge q at point f* to the dipole field at a 



\r - ro| J 

of a dipole moment K at point f . The monopole field 
is less singular at the location of the particle than the 
dipole field is. The former falls off as l/(|r — r~o|) while 
the latter falls off as l/(|f — ro|) 2 , in the limit r — > i%. 
Differently said, the field of the non-spinning particle is 
"less singular" in the limit r — > r o than the field of the 
spinning partilce is. 



V. FUTURE WORK 

The author and collaborators are currently working 
on a self-force calculation for a particle of mass m in the 
vicinity of a Schwarzschild black hole, as that calculation 
is described in [2(| • The ultimate goal is to generalize the 
calculation to the more interesting (for LISA data anal- 
ysis) case of a particle moving in the vicinity of a Kerr 
black hole. Knowledge of either the singular gravitational 
field given in Eq. (|79|) or of the direct contribution to the 
retarded gravitational field is necessary for those calcu- 
lations. 

The author is also working on calculating the part of 
the self-force that is due to the spin of a particle. The cal- 
culation of the singular gravitational field of a spinning 
particle will be presented in a future paper. Even though 
we intend to follow the general method described here, 
that calculation is more complicated than those described 
in this paper. At that time we will also present the spher- 
ical harmonic decomposition of the singular fields, which 
is necessaryfor the self-force calculation as was first de- 
scribed in [7| and outlined in Sec. IIIII in this paper. We 
will first carry out the spherical harmonic decomposition 
for the toy-problem of the scalar dipole field given in 
Eq. (f33|) to try and identify any issues that arise from 
the presence of intrinsic spin and then we will move on 
to performing the same calculation for the gravitational 
field of the spinning particle. 



Acknowledgments 

I wish to thank Steve Detweiler for his guidance in com- 
pleting this work and for his very useful suggestions on 
this manuscript. I also wish to thank Bernard Whiting 
for useful discussions. This work was supported in part 
by the Institute for Fundamental Theory at the Univer- 
sity of Florida and in part by NSF PHY-0200852 at the 
University of Wisconsin - Milwaukee. 



14 



LISA website available (http://lisa.jpl.nasa.gov/). 

B. S. DeWitt and R. W. Brehme, Annals of Physics 9, [18 

220 (1960). 

Y. Mino, M. Sasaki, and T. Tanaka, Phys. Rev. D 55, [19 
3457 (1997). 

T. C. Quinn, Phys. Rev. D 62, 064029 (2000). [20 
T. C. Quinn and R. M. Wald, Phys. Rev. D 56, 3381 
(1997). [21 
Y. Mino, H. Nakano, and M. Sasaki, Prog. Theor. Phys. 
108, 1039 (2003). [22 
L. Barack and A. Ori, Phys. Rev. D 61, 061502(R) [23 
(2000). 

L. Barack and A. Ori, Phys. Rev. D 67, 024029 (2003). [24 

L. Barack and A. Ori, Phys. Rev. D 66, 084022 (2002). [25 

L. M. Burko, Phys. Rev. Lett. 84, 4529 (2000). 

L. Barack and L. M. Burko, Phys. Rev. D 62, 084040 [26 

(2000). 

L. Barack, Phys. Rev. D 62, 084027 (2000). [27 
L. Barack, Phys. Rev. D 64, 084021 (2001). 
L. Barack, Y. Mino, H. Nakano, A. Ori, and M. Sasaki, [28 
Phys. Rev. Lett. 88, 091101 (2002). 

L. Barack and C. O. Lousto, PRD 66, 061502(R) (2002), 
gr-qc/0205043. 

S. Detweiler and B. F. Whiting, Phys. Rev. D 67, 024025 
(2003). 

E. Poisson, Living Rev. Relativity 7 (2004), gr- 



qc/0306052. 

S. Detweiler, E. Messaritaki, and B. F. Whiting, Phys. 
Rev. D 67, 104016 (2003). 

L. M. Diaz, E. Messaritaki, B. F. Whiting, and S. De- 
tweiler, Phys. Rev. D 70, 124018 (2004), gr-qc/0410011. 
T. S. Keidl, J. L. Friedman, and A. G. Wiseman (2006), 
gr-qc/0611072. 

K. S. Thorne and J. B. Hartle, Phys. Rev. D 31, 1815 
(1985). 

X. H. Zhang, Phys. Rev. D 34, 991 (1986). 
C. W. Misner, K. S. Thorne, and J. A. Wheeler, (Free- 
man, San Fransisco) (1973). 

S. L. Detweiler, Class. Quantum Grav. 22, S681 (2005). 
GRTENSOR software website available: (1992), 
http://grtensor.phy.queensu.ca/. 

J. Hadamard, (Yale University Press, New Haven) 

(1923). 

E. Messaritaki, Ph.D. Dissertation (2003), gr- 
qc/0308047. 

The differential equations become more complicated 
as the electromagnetic potentials and the gravitational 
fields are considered. For that reason, the differential 
equation is only shown for this simpler case of the scalar 
dipole field. It should be kept in mind that it has the same 
general form for the more complicated fields as well. 



